Chirality in Quantum Computation with Spin Cluster Qubits 
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We study corrections to the Heisenberg interaction between several lateral, single-electron quan- 
tum dots. We show, using exact diagonalization, that three-body chiral terms couple triangular 
configurations to external sources of flux rather strongly. The chiral corrections impact single qubit 
encodings utilizing loops of three or more Heisenberg coupled quantum dots. 
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Qubits based on electron spin in semiconductor struc- 
tures play a central role in solid state quantum informa- 
tion processing 0,0, LULU In particular, the extremely 
long p spin coherence times (T 2 ^> [is) coupled with the 
enormous lithographic and fabrication scalability advan- 
tages of existing semiconductor technology have led to 
a large number of proposed quantum computer architec- 
tures using electron spin qubits in semiconductor nanos- 
tructures. Most of these proposals [jj, L2, y, L4] utilize the 
Heisenberg exchange coupling ("the exchange gate" 0) 
between two electrons localized in individual quantum 
dots or in shallow donor states (e.g., P atoms in Silicon 
0) in order to carry out the required two-qubit opera- 
tions. The original proposals involve control, tuning, and 
manipulation of one-qubit and two-qubit gates through 
externally applied, local magnetic and electric field pulses 
respectively, but it was soon realized that controlling sin- 
gle localized electron spins through local magnetic field 
(and/or microwave) pulses is problematic [fjLijIa], per- 
haps even prohibitively so, and therefore a number of in- 
terestin g th eoretical proposals have recently been made 
La La La, Loi for using a cluster of spins ( "spin cluster" ) 
rather than a single spin, as the suitable qubit unit in 
solid state spin based quantum computing architectures. 
In particular, some elegant, recent theoretical literature 
EH EH EH 01 nas convincingly demonstrated as a 
matter of principle the non-essential and therefore dispos- 
able character of single-spin qubit operations altogether, 
showing that the Heisenberg exchange coupling between 
electron spins by itself can actually carry out quantum 
computation in spin-based architectures, provided each 
logical qubit is encoded in a number of physical spins 
(i.e., spin cluster encoded qubits). 

One important perceived advantage of such collective 
spin based logical qubits (avoiding single electron spin 
qubits altogether) is that the system can be placed in a 
decoherence free subsystem or subspace [3, EH (where 
the encoding scheme strongly suppresses decoherence of 
single logical qubits) making such encoded spin cluster 
qubits highly desirable from the perspective of scalable 
and coherent fault-tolerant quantum computation. Not 
surprisingly, therefore, the idea of using spin clusters 
rather than single spins as the basic building block of ex- 



change gate-based solid state quantum computation has 
attracted a great deal of recent theoretical attention. 

In this Letter we point out and theoretically formu- 
late an important conceptual, topological aspect of spin 
cluster qubits in solid state architectures. In particular, 
we show that the proposed solid state two (or three)- 
dimensional spin cluster qubits will necessarily generate 
topological chiral terms in the qubit Hamiltonian in the 
presence of external sources of flux, which will necessarily 
modify the Heisenberg exchange interaction bringing into 
question, in the process, the applicability of the exchange 
gate idea which has been the centerpiece of spin-based 
solid state quantum computation architectures. Thus the 
proposed spin cluster qubits involving looped arrays of lo- 
calized, tunnel coupled electron spins will not behave as 
originally envisaged in the exchange gate scenarios dis- 
cussed in many recent publications. Even the proposed 
encoded spin qubits, which specifically exclude all single 
qubit operations (and hence do not employ any applied 
external magnetic field), will have severe single-qubit de- 
coherence problems (thus undermining the decoherence 
free properties of the logical qubit) since any temporally 
fluctuating external sources of enclosed flux will cause 
dephasing through the chiral term discussed here. Con- 
tributions to the mean field component of an external 
magnetic field may arise, for example, from the nuclear 
dipolar field of the host lattice [5j . Alternatively, we find 
that the chiral term plays no role in clusters which do 
not contain three simultaneously, tunnel coupled sites. 
The chiral term may also be considered negligible when 
the area enclosed by three tunneling channels is small. 
Atomic clusters, with small lattice spacings ~ 1 A, re- 
quire large external mag netic fields, - 10 4 T, to see a 
sizable chiral contribution. We show that the large lat- 
tice spacing in looped quantum dot arrays enhances the 
role of three body chiral terms in the spin Hamiltonian. 

We study the low energy Hilbert space of N lateral 
quantum dots containing N electrons lying in the x — y 
plane by exact diagonalization of the following Hamilto- 
nian: 
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FIG. 1: Two possible looped configurations of lateral quan- 
tum dots in a magnetic field. The triangular configuration 
contains three tunneling channels with one loop formed by 
the vertices 123. The square configuration contains, in gen- 
eral, six tunneling channels with four three-site loops formed 
by the vertices 123, 234, 341, and 412. 

where we define the canonical momentum, P = p + -A, 
and the Zeeman term Hz — g*fisS ■ B. In GaAs we 
have an effective mass m* = 0.067m e , dielectric con- 
stant e = 12.4, and (/-factor g* = —0.44. We work in 
the symmetric gauge with magnetic field B — Bz. The 
field couples directly to the total spin, S, through the 
Zeeman term. The single particle potential confines the 
electrons to lie at the vertices of an equilateral triangle 
and square for N = 3 and 4, respectively, as shown in 
Fig. H In diagonalizing the N = 3 system we choose: 
V"(rj) = (m*Wg/2) minj{|ri — Rj| 2 }, where a parabolic 
confinement parameter, hu>o = 6 meV, effectively local- 
izes the electrons at the sites Rj separated by 40 nm. 

Prior to exact diagonalization of Eq. we seek a per- 
turbative expansion in terms of on-site spin operators. 
We consider the single band, tight binding limit, a good 
approximation in the limit that the excited states of the 
quantum dot lie higher in energy than the lowest spin 
split states. We also take the on-site Coulomb interaction 
to be much larger than the tunneling energy. (We verify 
numerically that, for the N = 3 system, the following ap- 
proximations are self consistent only at low B, oj c < Wo/3, 
where uj c = eB/m*c is the cyclotron frequency.) As a re- 
sult we obtain the Hubbard Hamiltonian: 

Hn = -^2 tij4a c jc* + ( Un n n il + 9*^bB ■ Si) , (2) 

where c ia creates a fermion at the site i with spin a €t! 
and rti a — c\ a c ia . We have defined the spin operators 

S< = ? c L' T aa' c ia' wnere c are the Pauli matrices. As 
a consequence of the magnetic field the tunneling coef- 
ficients are complex: £y = |ty| exp(27ri$y /$o), where 
<J>o = hc/e. The magnetic vector potential generates the 
Peierls phase: $^ = J? A • dr. The integral is taken 
along a path connecting the sites i and j. Working in 
the small \tij\/U limit we confine our attention to single 
occupancy states. Here we expand Hn by applying a 
unitary transformation exp(iK)Hn exp(— iK), where K 
is an operator changing the number of doubly occupied 



states. Up to third order in this expansion we find |l7lll8| : 

i i,3 

+ J2 ^S r (S J xS,)+^(^], (3) 

where t ~ The second term is the usual Heisen- 

berg interaction, where Jij = 2\tg\ 2 /U, while the third 
term is a three site sum over chiral 19] terms Xijk = 
MyfeSi • (Sj x Sfe) around loops, A. In H ff we have ex- 
cluded fourth order terms of the form (Sj X Sj) • (Sfe x S;) 
and (Si ■ Sj)(Sfe • S;). It was shown that the later plays 
a role at B = in a four spin, tetrahedral configura- 
tion [2(|. The coefficients in the chiral term: /i.yfe = 
(24/C/ 2 )|^||^fe||ife l |sin(27r$yfe/$o), depend on the flux 
enclosed by the three site loop, <&ijk- As a result, the 
chiral term generates an energy splitting between third 
order, virtual tunneling processes which run along and 
counter to the applied vector potential. The phase, 
27r$yfe/$o, is the Aharonov-Bohm phase generated by 
the virtual current moving around the flux enclosed by 
the three site loop. The chiral term is Hermitian but 
breaks time reversal symmetry and vanishes on bipartite 
lattices as a result of particle hole symmetry 0] . 

We may now ask whether or not three coupled quan- 
tum dots containing three electrons may support a noise- 
less subsystem in the presence of a fluctuating, perpen- 
dicular magnetic field |l5| . The set of states comprising a 
noiseless subsystem remain invariant under the applica- 
tion of a suitable noise operator. We may, for example, 
choose the first term in Eq. @ as a noise source. We 
then find the simplest example of a quantum dot, noise- 
less subsystem in the 5=1/2 sector of the N — 3 system. 
This encoding makes use of a fourfold, B = degener- 
acy found in this sector to protect quantum information 
stored in the qubit defined by |A)a, where A = or 1. 
The subscript denotes the number of quantum dots and 
electrons used to define the qubit. The four states may 
be written: 

ia> 3 ® i-i) = ^(mt)+^ +1 im)+^- A itu» 
ia> 3 ® |i) = -i=(im)+^ A+1 inT) + ^ 2 - A uTT))(4) 

where to = exp(27ri/3). In the presence of a fluctuat- 
ing Zeeman energy the first term in the tensor product 
preserves the quantum information stored in the quan- 
tum number A while the z-component of spin, the second 
term in the tensor product, may fluctuate. Application 
of H c s up to second order shows that the magnetic field 
only weakly affects |A)3 through the Heisenberg term. In 
fact several proposals 0, El G3 El E3] suggest use of 
the Heisenberg term, with anisotropic couplings Jy , to 
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implement Pauli gates on encoded qubits. In Heisenberg 
gating schemes the expectation value of the individual 
Heisenberg terms, and hence gates formed from them, 
remain insensitive to fluctuations in the Zeeman term 
because [Si ■ Sj,^2 i Si] = for all i and j. 

Up to second order, H c g also allows an encoding 
against Zeeman-like or collective noise. Here collective 
noise implies symmetry among all spins when coupling 
to a bath. However, the third order chiral term acts as 
a non-collective system-bath interaction between the en- 
coded qubit and a potentially noisy source of enclosed 
flux. The chiral and Zeeman terms remove all degen- 
eracies required to establish a qubit encoding immune to 
fluctuations in the perpendicular magnetic field. Explic- 
i%: Xi23|A) 3 ®|±|)=Mi23(2A-l)(N/3/4)|A) 3 ®|±i), 
where /ii2 3 = (12t J/U)sin (27r$i23/$o) hi the case Jy = 
J and \tij\ = t for all i and j. Therefore, X123 yields 
an effective Zeeman splitting between the encoded basis 
states of the three spin qubit. The size of the splitting 
depends on 6J(J/t) which, under reasonable conditions, 
cannot be neglected as long as the second order term, 
JijSi • Sj, remains large. Gates constructed from the ex- 
change interaction must operate on time scales shorter 
than typical spin relaxation times supplying, in turn, a 
necessary minimum to the exchange interaction. 

To test the accuracy of H e g we diagonalize Eq. (JTJ 
with N = 3. We construct the matrix representing H 
in the Fock-Darwin basis centered in the triangle 
formed by the three parabola centers defined by V(r). 
We find it necessary to include up to ~ 10 5 origin cen- 
tered, Fock-Darwin basis states with z-component of an- 
gular momentum less than twelve to obtain convergence. 
We use a modified Lanczos routine to obtain the ground 
and excited states. We focus on the three lowest energy 
states in the absence of the Zeeman term. Fig. [5] shows 
the energy of the lowest states obtained from exact diag- 
onalization of H in the S z — 1/2 sector as a function of 
magnetic field. The ground state energy is set to zero. 
The two lowest energy states have total spin S = 1/2. 
They are degenerate at B = as expected from the re- 
flection symmetry of the triangular confining potential. 
The next highest state has S = 3/2 which corresponds 
to 6t 2 /U w 0.125 meV. Above this state we find (not 
shown) the higher excited states to lie near 3.8 meV. 

As we increase magnetic field the magnetic vector po- 
tential breaks the reflection symmetry of the confining 
potential leading to a splitting between the two lowest 
states. The splitting is linear in B, for small B, which 
suggests that the splitting is due to Xi23- The chiral term 
annihilates 5 = 3/2 states. We therefore expect that the 
energy of the S = 3/2 state, E s=3 / 2 , does not change 
with magnetic field at low fields (while Eg =3 / 2 — -E|o> 3 
should increase linearly). Here the relevant length scale is 
a modified magnetic length, a — y/hc/eB(l+4uj 2 /Lj 2 )^i . 
As a result the magnetic field increases the Coulomb en- 
ergy (~ e 2 /ea) only above B ~ 0.7 T. At fields above 
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FIG. 2: The energy of the lowest states of three electrons 
confined in three adjacent, parabolic quantum dots obtained 
from exact diagonalization of Eq. fl) . The dot centers form 
an equilateral triangle in the x — y plane with 40 nm side 
lengths. The energies are plotted as a function of perpendic- 
ular magnetic field. The total z-component of spin is taken 
to be 1/2. The energy of the state |0}3 in Eq. J2J is set to 
zero. The solid and dotted lines show the energy of the states 
with total spin, 5 = 1/2 and S = 3/2, respectively. The 
dashed line plots the Zeeman splitting in GaAs, w 0.025B[T] 
meV, for comparison. The energy splitting between the three- 
electron encoded qubit states 1 0) 3 and 1 1 ) 3 is larger than the 
Zeeman splitting between the spin up and down states of a 
single electron. 



B ~ 0.7 T the energies slope down indicating an even- 
tual sign change in J, as for double quantum dots |23| . 
Here the chiral contribution starts to become suppressed, 
along with the exchange interaction. At larger fields the 
long range part of the Coulomb interaction becomes rel- 
evant. We must then keep extended terms in of the 



form V Y^, 



n-i+i.a' ■ We note that inclusion of 



these terms into Hr merely renormalizes the on-site con- 
tact term U for U ^> V 24]. For large V, the extended 
Hubbard term favors double occupancy of the dots and 
eventually leads to a sign change in J 23] . This suggests 
that H e g is qualitatively accurate below B ~ 0.7 T. 

The slope of the energy splitting between the two low- 
est states in Fig. [5] allows us to estimate t/U for this 
system. Using the area of the equilateral triangle defined 
by the dot centers we obtain t/U « 0.09 which shows that 
our expansion in t/U is consistent at low fields |2j|. For 
N = 3, only odd powers of allow linear magnetic field 
dependence in the splitting, showing that the magnetic 
field dependence captured by the chiral term is accurate 
up to <d 

The dashed line in Fig. [2] shows, for comparison, the 
Zeeman splitting between the spin states of a single elec- 
tron. The Zeeman splitting is smaller than the energy 
difference £ , |i) 3 — E\o) 3 suggesting that, in the chosen 
parameter regime, quantum information encoded in the 
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"noiseless" spin states of a three electron-triangular set of 
quantum dots is more sensitive to a perpendicular mag- 
netic field than a single electron. 

We now study magnetic field effects on a decoherence 
free subspace (DFS) [l4[ formed from quantum dot loops. 
The lowest number of physical spins supporting a DFS 
is four. For simplicity, we begin with four quantum dots 
containing four electrons lying at the vertices of a square 
with equal tunneling |iy| = t, including diagonal terms 
shown schematically in Fig. ^ In this case we find a DFS 
among two S = states corresponding to A = and 1: 



|A)4 = imi} + iun)+^ A+i inu) 
mt)+^- A i mi)+^ 2 - A i tut). 



,A+1| 



(5) 



These states have equal Zeeman energies and, up to sec- 
ond order in Eq. 0, show no direct magnetic field depen- 
dence in their energy spectra. As for the third order term 
the spin Hamiltonian must respect the exchange sym- 
metry inherent in the lattice. Therefore, the sum over 
asymmetric chiral terms must vanish with equal tunnel- 
ing among all sites. In the |A)4 basis we find: 

^ X«fc ® ij|A>4 = -2~(2A- 1) ^2 M»ifcey«|A>4, (6) 
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where tijki is the four component Levi-Civita symbol and 
Ii is the identity operator. The sums run over three site 
loops excluding I = i,j or k. The sum vanishes with 
tunneling |ty| = t for all % and j. However, if we impose 
some asymmetry in tunneling, as is done when applying 
a Pauli gate to the encoded qubit, this sum is not zero. 
Consider a simple case: |i 3 i| = |t 2 3| = ^34] = t(l + 
S), where S is a number and all other |fy| = t. The 
sum over chiral terms then induces an energy splitting 
~ 24irV3tJ5AB z /(U$ ) between the states with A = 
and 1 for <f>yfe/$o "C 1. Here A is the area of the triangle 
defined by the vertices 123. As in the N — 3 system, the 
energy splitting acts as an effective Zeeman splitting on 
the N = 4 encoded qubit with the exception that here 
the parameters A and S may be included in an effective 
^-factor of the encoded two-level system. 

As mentioned earlier, Pauli gating sequences may be 
applied to encoded, multi-spin qubits by tuning the 
Heisenberg couplings, and therefore fy, in an asymmet- 
ric fashion. When applied to a DFS a Pauli gate com- 
posed of Heisenberg terms must, by construction, involve 
a spin specific asymmetry with, for example, large 6. As 
shown above, the N = 4 encoded qubit will, during a 
gate pulse, be sensitive to sources of enclosed flux be- 
cause the scalars forming the Heisenberg interaction do 
not commute with the individual chiral terms where, for 
example: [Si ■ S3, X123] 7^ 0, with i = 1 or 2. 

We have shown that, in multiple quantum dot archi- 
tectures containing tunnel coupled loops, fluctuations in 
enclosed flux provide a potential source of phase flip er- 
ror in noiseless subsystems and subspaces through chiral 



currents generated by virtual hopping processes around 
three site loops. In systems for which fluctuations in the 
phase of tij can be ignored we have shown that a well 
controlled source of flux may also be a good candidate 
for implementing a logical Pauli Z gate on a three or 
four spin encoded qubit. Quantum algorithms using only 
exchange based quantum gates require several, accurate 
applications of the exchange Hamiltonian Q. The chi- 
ral term may then supplement the exchange interaction 
and therefore potentially reduce the overhead in exchange 
based algorithms. However, the phase in tij is in gen- 
eral a Berry's phase |2(| arising from symmetry breaking 
terms in the original Hamiltonian. Spin orbit coupling 
can contribute to Berry's phase effects in triangular lat- 
tices [22! and may therefore play a role in effective spin 
Hamiltonians modeling tunnel coupled quantum dots. 

This work is supported by ARDA and NSA-LPS. 
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